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I. INTRODUCTION 

Over the past few years, the holographic gauge-gravity duahty has gone through a series of 
transformations, moving further away from its most comfortable arena, that of theories with 
a large amount of supersymmetry. Top-down studies, based on string/M-theory constructions, 
have been met by a number of bottom-up approaches, with applications ranging from QCD-like 
theories to condensed matter systems. For example, holographic techniques have been applied 
to probing the transport properties of the strongly coupled QCD quark gluon plasma (QGP), 
with the notable order of magnitude agreement of the shear viscosity to entropy density ratio ^ 
between RHIC (and now LHC) data and the holographic result first found in [Ij^. 

In the context of applications to condensed matter physics, recently much focus has been on 
non-relativistic systems undergoing quantum critical behavior, which often exhibit an anisotropic 
scaling symmetry of the type 

t^X^t, x->Ax, (1) 

with the dynamical exponent z generically not equal to one. On the gravity side of the duality, 
this scaling can be realized as an isometry of the metric 

ds^ = r^-dt' + r^dx^ + '^ (2) 

as long as one also sends r — )• X~^r. For the case of pure AdS (in the Poincare patch), the 
dynamical critical exponent is simply z = 1, describing relativistic scaling. Studies of holographic 
duals exhibiting non-relativistic symmetries were initiated in [6l [7] , with the first examples of 
Lifshitz scaling found in [8j. There exist several string theory constructions of non-relativistic 
holographic duals, including many which exhibit Lifshitz scaling [9ll24j. There is by now an 
extensive literature on the subject, and we refer the reader to the review articles |25f[28j for a 
more detailed treatment of holographic applications to condensed matter physics. 

Our interest in systems exhibiting the anisotropic scaling symmetry ([T]) stems from a rather 
different reason. We would like to understand what are the ingredients needed to generate 



^ For applications to rj/s see for example the review articles [IHl], and for an exhaustive review of the applications 
to the QCD plasma see and references therein. 
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temperature dependence for the shear viscosity to entropy ratio r//s of strongly coupled field 
theories, within the framework of the holographic gauge/gravity duality. It is by now well- 
known that in all gauge theories with Einstein gravity duals^ one finds the universal result 
rj/s = l/47r \34:\ I35j. Moreover, as argued originally via the membrane paradigm in [36j, t//s is 
given strictly by horizon quantities, a result which holds even in theories with higher derivatives 
|37H39j . Although r]/s does not run in any Wilsonian sense^ (the associated radial flow on the 
gravitational side is trivial, explaining why it can be extracted from horizon data) it can still 
undergo temperature flow. As a concrete example of a setting where such a flow is generated, we 
will consider a system which exhibits non-relativistic scaling at low temperatures, and standard 
relativistic scaling in the high temperature regime. We should note that having a handle on the 
behavior of 7//s as a function of temperature is relevant for the physics of the strongly coupled 
QGP, and in particular for better understanding the elliptic flow measurements from heavy ion 
collisions at RHIC and LHC (see e.g. milS]). 

Our starting point will be Einstein-Maxwell-dilaton gravity in five dimensions. Thanks to 
the presence of the dilatonic scalar field, finite temperature solutions in this theory (with a finite 
charge density) exhibit Lifshitz-like scaling near the horizon, while asymptoting to AdS at large 
distances [36l - [5l] . Of particular interest to us, however, is the temperature dependence of the 
solutions. In fact, the Lifshitz-like scaling regime appears in the near-horizon geometry when 
the temperature is much smaller than the chemical potential set by the finite charge density, 
T jjL [i.e. in the IR). On the other hand when T ^ fi (i.e. in the UV), the near-horizon 
solution is replaced by an anti-de Sitter black brane. The interpolation between the low-T and 
high-T regimes has been shown to be smooth - no discontinuous phase transition occurs as T is 
varied [50\ [5T] . An interesting temperature flow for the viscosity to entropy ratio is generated 
by adding appropriate higher derivative corrections to the theory - without which there would 
be no deviation from the universal r^/s = \/Att result. 

^ An exception is the case of anisotropic fluids. A way to obtain a deviation from universality by breaking 
the rotational symmetry spontaneously was shown in [291 I30| . Moreover, the resulting non-universal shear 
viscosity to entropy ratio has been shown to exhibit an interesting temperature flow, without the need for 
higher derivative corrections. This can be seen in '3D -32j in the context of p-wave superfluids, as well as in |33) . 
in the anisotropic axion-dilaton-gravity background of |23l I24| . 
For recent attempts at refining notions of holographic RG flow, see [40H43| . 
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Once such higher derivative terms are turned on, we will find that the behavior of r//s in the 
T <^ fx regime is sensitive to the dynamical exponent z, while for T ^ /i it is entirely independent 
of it - knowledge of z has been washed out by the large temperature in the system. Moreover, 
r]/s will interpolate smoothly between its IR and UV values as the system is heated up. While 
this is not Wilsonian RG flow, it does offer an example of a setting in which the structure of 
r]/s in the IR is very different from that in the UV, and in which UV data alone is not enough 
to specify the behavior of rj/s. Another example of such a stark difference between UV and IR 
physics - and of an analogous temperature flow - was seen in |55j , where the system underwent 
a second order (superfluid) phase transition^. Just like in the model of [55], here the non-trivial 
behavior for the viscosity to entropy ratio can be traced to the presence of a matter field, and 
in particular to its coupling to the higher derivative terms. In fact, as we will see, it will be the 
dependence oi r]/s on the (horizon value of the) scalar which will lead to additional temperature 
dependence, and a significantly more intricate behavior in the IR. Although a dependence of 
r]/s on T/fi was already seen in theories with higher derivative corrections at finite chemical 
potential [37^ 156) . its structure was not as rich as the one observed in this setup, where it is 
sensitive to the dilatonic couplings and the different scalings in the low-T and high-T regimes. 

To summarize, our main goal in this paper was to gain insight into what features are needed 
to construct a system in which transport coefficients such as rj/s undergo temperature flow. 
The model we have chosen here provides a concrete example - in no way unique - of such a 
setting, where the evolution of a scalar field translates into a flow for tj/s. Interestingly, in this 
model a particularly non-trivial temperature dependence arises when a shift symmetry of the 
two- derivative action is explicitly broken by the higher-derivative terms - linking a much richer 
behavior of the transport coefficients to the presence of broken symmetries. In fact, breaking 
of this shift symmetry will be crucial for achieving a particularly significant deviation from the 
behavior of the standard Reissner-Nordstrom anti de-Sitter (RN-AdS) black brane background 
studied in [371 [56]. Finally, we should mention that, although understanding the intermediate 

■* In the construction of [55] the higher derivative terms were engineered to be non-vanishing only below the 
temperature Tc associated with the phase transition. As a result, in 'SF one finds the universal result rj/s — l/4n 
above Tc, while a deviation from it below Tc, with an associated running with temperature. 
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temperature regime (when T ~ /i) in our model would be interesting, such an analysis is beyond 
the scope of this paper, and we have restricted our attention to the endpoints of the flow. 

The organization of the rest of the paper is as follows. In Section |II] we introduce the 
Einstein-Maxwell-dilaton theory we will focus on, and discuss relevant properties of its black 



brane solutions. We present the details of the calculation of rj/ s in Section III We also include 
an expression for yy/s in terms of a generic near-horizon black brane expansion, which will 



facilitate the discussion of temperature dependence. Section IV is dedicated to the flow of ij/s 
as a function of temperature. We discuss our main results in Section IVl 



II. THE SETUP 



We take as a starting point the following Einstein-Maxwell-dilaton action, 

So = —V / d'xV^iR -2A-2 {V<pf - e^'^'f'F^) , (3) 
lb7rG5 J 

where (p is a real scalar and F = dA is an abelian field strength. Thanks to the presence of the 
U{1) gauge symmetry in the bulk, the dual field theory has a conserved particle number, with 
an associated chemical potential /x. The action is also invariant under the shift symmetry 

(j)^(l) + 6, e-^^A , (4) 

reflecting the fact that the normalization of the gauge coupling can be absorbed into the defini- 
tion of the gauge field (and in particular, in the definition of the charge). 

Since we are interested in higher-derivative modifications to the shear viscosity to entropy 
density ratio, we will add to the leading order action Sq curvature corrections of the form 

Si = I d^'x^g (/3i e^i'^iJ^.p.i?'^'^^'^ + (32 e^^'^iJ^.p.F^'^F'"^) , (5) 

with 7i,72 arbitrary constants. We will take the couplings (3i and /32 of the higher derivative 
terms to be perturbatively small, as would be the case in a derivative expansion within a genuine 
string theory reduction. Note that generic scalar couplings 7i 7^ explicitly break the shift 
symmetry Q of the two-derivative action. However, for the special case of 71 = and 72 = 2a 
the symmetry is preserved. Whether the shift symmetry is broken or not will play an interesting 
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role in the discussion of t//s - and in particular of its temperature dependence - as we will see 
in Section IV Finally, although ([s]) is clearly not the most general form for the action at the 



four-derivative level, the terms which we have included are the only ones that will contribute to 
r//s, as we have checked explicitly^ in accord with the general arguments put forth in [37J. 

Before turning to the discussion of higher derivative corrections to r]/s, we would like to 
highlight some of the properties of the solutions to the leading order two-derivative action 
Sq. We will restrict our attention to electrically charged black brane solutions to ([3]), whose 
thermodynamics has been investigated recently in a number of studies |46l - [5^ , aimed at applying 
holographic methods to strongly coupled condensed matter systems.^ The finite temperature 
solution to ([3]) which is perhaps best known is the uncharged AdS black brane, 

ds^ = -f{r)dt'+'^^dx' + ^^, m='^J^-(-y) , (6) 



0(r) = 00, A = Aodt, (7) 

in which the gauge field and the dilatonic scalar field are both constant. Here the radius L 
is determined by the size of the cosmological constant, A = —Q/L'^. For this solution the 
temperature and entropy density are given by 

T=— s= = (8) 

However, what we are interested in are solutions with a finite charge density^ and therefore a 
non-trivial gauge potential. To this end, we note that the action ([s]) also admits a black brane 
solution with a metric exhibiting the anisotropic Lifshitz scaling ([T]), a non-trivial gauge field 
and a logarithmically running dilaton. The scaling symmetry in this case is not exact, but is 
broken by the running of the dilaton - as a result, the solution has been termed Lifshitz-like 
|47j . Taking the following ansatz for the fields, 

ds"^ = -a{rfdt^ + b{rfdx^ + c{rfdr^, 

</. = (/)(r), A = A{r)dt, (9) 



^ We note, however, that we have not considered the addition of terms involving derivatives of the curvature, 
such as e^^'^V^i? or 6^**^X7" V*" -Rat (or terms related to these by integration by parts). Their inclusion is more 
subtle, since the application of the Wald entropy formula in this case becomes slightly more involved |57) . 

^ See also [58] for a discussion on the thermodynamics of Lifshitz geometries including higher derivatives. 
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the Lifshitz-like black brane solution to the action ([s]) is given by |46fl48t [53] : 

The parameters specifying the solution are determined by the charge Q of the U{1) gauge field 
and the dynamical exponent z, which is itself fixed by the value of the dilatonic coupling a : 

12 ' 6L4(z-l)' a2 ^ ^ 

Note that with our notation the dimension of the charge is [Q] = L? . Also, the temperature of 
this scaling solution is 

with the appropriate units of energy, [T] = 1/L. Noting that the entropy density still behaves 
as s ~ rj^ and that the chemical potential also has units of energy [/i] = 1/L, simple dimensional 
analysis tells us that the entropy must scale as s ~ . 

In the Lifshitz-like geometry the gauge field and the e^^'^F^ term in the action diverge in the 



large r region. The divergence can be cured, however, by embedding the scaling solution (10) 
into an asymptotically AdS geometry, with the embedding providing a relativistic UV completion 



of the model. The new geometry can be obtained by adding to the scaling solution (10) a 
perturbation that is irrelevant close to the horizon, but which becomes increasingly important 
for larger values of r. In the new asymptotic region both the scalar and the gauge field approach 
constant values - their divergence effectively cut-off by the scale at which the theory becomes 
AdS-like. It is the chemical potential /i which plays the role of the cut-off, setting the scale for 
the cross-over between the two distinct behaviors of the geometry. Finite-temperature brane 
solutions interpolating between a near-horizon Lifshitz-like regime and an asymptotically AdS 
region - realizing the embedding described above - have been constructed explicitly in |47H51j . 

It is important to emphasize that in the dilatonic system described by ([s]) the nature of the 
solutions is sensitive to the hierarchy between the two scales in the theory, the temperature T 
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and the chemical potential fi. At a fixed temperature T <^ ^, the near-horizon geometry r ^ 



is that of the Lifshitz-like brane (10), while the large r behavior of the metric and gauge field is 



described by the same asymptotics as those of a planar charged AdS black hole, namely, 

*('■) = 00 + ^ + ..■ , (13) 

where the ellipses denote terms of higher order in 1 /r^ . The physical chemical potential /i can 
be read off as usual from the non-normalizable mode of the gauge field, and we have allowed 
for an additional constant e"""^" for reasons which will become clear shortly''. Away from the 
horizon and the boundary, the full geometry does not admit a simple analytical form and one 
must resort to studying it numerically (we refer the reader to [371 SSHSI] for a discussion of the 
numerics) . 

On the other hand, when the temperature is the dominant scale in the system, T ^ ^, 
the geometry is AdS-like everywhere. In fact, when T is so high that the charge is essentially 
negligible, we can think of the uncharged AdS brane ^ as describing the entire geometry, from 
the horizon to the boundary. With the exception of this very large temperature regime, the full 
solution for T > fi should be thought of as a perturbation of whose behavior at intermediate 
distances is again only known numerically. Solutions in the intermediate temperature region have 
been analyzed numerically in |50tl51j. where it was shown that the system evolves smoothly from 
the low-T to the high-T regime. 

In summary, in the construction sketched above we see two types of 'flows.' One is the stan- 
dard radial flow between a near-horizon geometry (Lifshitz-like if T ^ ^) and AdS asymptotics. 
The other flow occurs as a function of temperature. By heating up the system, the near-horizon 
Lifshitz brane which is present at T ^ is replaced by an AdS brane, with the chemical poten- 
tial controlling the onset of the transition. Once the chemical potential is negligible compared 
to the temperature of the system, T ^ /x, the AdS geometry is the one that dominates. As we 



We could have absorbed it into the definition of charge. 
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already mentioned above, the interpolation between the two opposite regimes is smooth |501l51j. 
It is the latter temperature flow that we will be mostly interested in, since it will provide us 
with precisely the type of setup we need to probe the behavior of r//s. Also, while it would be 
interesting to study the system at arbitrary temperatures, here we will focus only on the two 
limiting cases T <^ fi and T ^ fj,, which will be sufficient to illustrate our main points. 



III. SHEAR VISCOSITY TO ENTROPY RATIO 

We are now ready to compute the viscosity to entropy ratio^ for the full theory 

^ f d^x^/^ \r-2A- 2{V(t^f - e^^'^F^ 
vrGs J L 



under the assumption that the higher derivative couplings /3i and f32 are perturbatively small. 
Throughout the paper we will be working to linear order in the couplings /3j. For reasons 
which will become clear shortly, for the r]/s calculation it is particularly convenient to make the 



coordinate change u = rf^/r and rewrite the black brane solution (10) in the form 



where for the Lifshitz-like brane one has: 



A = h{u) dt, e2"<^(") = , (15) 



f{u) = 1 - u^^^^y' , h{u) = ^^^^ n-(^+3)/V(n) . 



(16) 



The boundary and the horizon are now located at n = and n = 1, respectively. We note that 
the standard AdS brane solution ([g]) also fits into this form, with z = 1 and the scalar field and 
gauge field set to fixed values. 

Although the metric ( |15[ ) is a solution to the two-derivative theory ^ only, we emphasize 
that for the computation oi-q/s knowledge of the back-reaction due to ^ is not needed. This is 



Dilatonic couplings to higher derivative terms in the context of rj/s have also been explored in |59| . We should 
also point out that Gauss-Bonnet corrections to rj/s in the model described by our two-derivative action Q 
have been studied in (48j where, however, the higher derivative terms were not coupled to the scalar. 
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a simple consequence of the universality of the shear viscosity to entropy ratio in asymptotically 
AdS backgrounds. To see why this is so, recall that for Einstein gravity r]/s = l/47r, indepen- 
dently of the form of the solution. The terms that would lead to a deviation from this universal 
result can only come from the higher-derivative corrections in the action, which are already of 
order 0{f3i). Thus, working to first order in the perturbative parameters /3j, one only needs the 
background solution - the back-reacted solution will yield a correction which is second order in 
perturbation theory. In particular, this means that when computing ij/s we don't have to take 
into account the renormalization of the dynamical exponent z found in [60] . 

In the framework of the gauge/gravity duality, the shear viscosity can be computed in a 
number of ways. Here we will make use of Kubo's formula, which relates r/ to the low frequency 
and zero momentum limits of the retarded Green's function of the CFT stress tensor: 

ri = -lim-ImG^ {uj,k = 0). (17) 

Let's sketch the main steps that go into finding r//s, following closely the discussion^ in [371 164| . 
The retarded Green's function can be extracted from the effective action Seff for the shear 



metric perturbation hx^{t,u) = f p^0fc(u)e Expanding (14) to quadratic order in 

the fluctuations (f)kiu) gives the by- now standard form for the effective action |64] . 



5. 



// -fl^du [A{u) <p'lct)_k + B{u) ^'k<j)'_k + C{u) + 

+ D{u) (^k^^k + E{u) + F{u) 0'^</.'_,] + Sgh , (18) 

where Sgh denotes the (generalized) Gibbons-Hawking boundary term^*^ and the coefficients 
A{u), B{u), . . . , F(u) encode information about the background geometry. After a series of 
manipulations the effective action can be shown to reduce to the boundary term 



Seff = J J^kllJo ' (19) 

with the flux term J^k given by 

2J-fc = Uk(l>-k + (C - A')^k(l)~k + ^0fc</''-fc + . • • , (20) 



Alternatively, rj/s can be extracted from the effective graviton coupling [61tl63l 



9 

For the explicit expression for the Gibbons-Hawking boundary terms we refer the reader to [HJ. For some of 
the subtleties involved in constructing boundary terms and counterterms in theories with curvature-squared 
corrections in the presence of a chemical potential see e.g. |65] . 
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where we have omitted a contribution coming from the generahzed Gibbons-Hawking boundary 
term^^, and Ilk(,u) is the radial momentum conjugate to (j)k- 

U,iu) ^ ^ = ((B - ^ - lF')<p',iu) - {E<p'au)y^ . (21) 

The retarded Green's function can be shown to be related to the flux through: 

^^y^^y = - ^1^0 Mu)7-kin) ■ ^^^^ 



Moreover, the only term in the flux that is relevant for the computation (17) of the viscosity 



can be shown to be the first, 2Tk = 4'-k + • • • • This, combined with (22), implies that the 



Kubo relation (17) reduces to the simple form 



ri = hm ."-'^=°^;\ , (23) 
relating the shear viscosity to the boundary behavior of 11^ and <j)k, in the low- frequency limit 

m- 

However, one can do even better. In terms of the radial momentum, the equation of motion 
for the fluctuations (pk can be recast in the following way, 

duUk{u) = M{u)(l)k{u), (24) 

where M(u) plays the role of an 'effective mass' for the fluctuations, and is given by: 

General arguments [37] show that the effective mass must vanish in the low-frequency limit"*^^, 

M{u) = Oioj^) . (26) 



We have verified this explicitly for our four-derivative action (14) by evaluating M{u) using 



the generic ansatz confirming that in the hydrodynamic regime the momentum Ilk{u) is 
These terms will not be needed here. 

Here we briefly outline the argument of [37]. Consider the case of a constant (/I>(it), for which the momentum II 



vanishes automatically, as clear from ([21|. Expanding the mass term as a series in the frequency, one can then 
show that to satisfy the equation of motion we must have at least M{u) = C'(cj). Time reversal invariance then 
imposes that M{u) = 0{lo'^), showing that the effective mass vanishes in the hydrodynamic limit. 
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independent of the radial direction u. Similar arguments [37j show that, to leading order in the 
low-frequency limit, uj (j){u) is also independent of the radial position^^. Thus, we are free to 



evaluate (23) at any value of u, and in particular at the horizon radius, which simplifies greatly 



the computation of ij/s in many instances. What we have seen, then, is that the shear viscosity 
is intrinsically a near-horizon quantity [3(ij, even in theories with higher derivatives. 



Putting all the ingredients together, the final expression for r] given a metric of the form (15) 



and the effective action (18) can be written [37] in the compact form 

1 



rj 



87rG5 

where Uh denotes the horizon radius, and 



K2{uh) + KiiUh) 



K2{u) 
K4(n) 



/ guu[U) 

9tt{u) 
E{u) 



A{u) - B{u) + 



F'{u) 



9uu 

Qtt 



For our black brane metric (|15|) we arrive at: 

2(32; 



1 + 



Pi 



I) 



levrGsL^ 1 AL2 

To find the entropy we make use of the standard Wald entropy formula [66! 

5C 



-27r / d^x^ 



-h 



(27) 

(28) 
(29) 

(30) 

(31) 



with h the induced metric on the horizon cross section S, and e^u the binormal to S. For our 



action and for the metric (15) we find that the entropy density s = S/oj^ is 

2 



1 + 



AL2 



(2(3z - 2)f{uh) - Anuh))e^^'f"^f3i + - (h'^f e^^*^/32 \, (32) 



where denotes the (infinite) black brane area. Finally, the expression for the ratio of shear 
viscosity to entropy density is 

2 J, 



s 



1 

An 



(33) 



At two-derivative level this follows from the fact that 11 is held fixed in the low frequency limit and so uj(f>'{u) oc 
cjn(u) vanishes. The fact that it still holds in a perturbative expansion in the higher-derivative couplings simply 
follows from the two-derivative result |37| . 
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parametrized in terms of the values of f{u), h[u) and (t){u) at the horizon, as well as the couplings 
of the higher-derivative corrections in the theory. 

A. Computing 77/5 for a generic near-horizon solution 

Given that the shear viscosity is a horizon quantity, as we have discussed above, we would like 
to close this section by showing how r}/s depends on the parameters which specify the horizon 
data of a black brane solution^^. This will prove useful in the next section, once we discuss the 
behavior of rj/s as one varies the temperature of the system. With this motivation in mind, 
we start by writing down a generic near-horizon expansion for our scaling solution, assuming a 
first order zero in gu and a corresponding first order pole in g^r, and imposing regularity for the 
gauge field at the horizon. Under these assumptions, the near horizon expansion is then of the 
form 

a{uf = ao(l -u) + ai(l - u)^ + 02(1 - uf + ... , 
b{uf = bo, 

C{uf = Co(l - U)"^ -I- Cl -I- C2(l -U) + ... , 

A{u) = Ao{l-u) + Ai{l-uf + ..., 

<P{u) = (f>h + M'i^ - u) + 02(1 - uf + ... . (34) 



Feeding this type of expansion into the prescription (27) we find that the shear viscosity is 

_ V ft (aoco-aoci + 3coai + 2aoco0i7i) 2fl -fi6^\ foK\ 
levrGs V ooCq / 

unaffected by the R^^paF^'^ PP" term in the action, in agreement with |37t I56| . The entropy 
density can be shown to be given by: 

^ ^ / _ (3coai - aoci) ^2^^^,.^, _ 2A| ^2 (gg) 
4G5 V ao% aoco J 

Combining these two expressions, we find that the viscosity to entropy density ratio reduces to: 

1 = ^(1- (1 + 27101) ^2^^^,.^. + ^L2/32e^^^'^') . (37) 
s 47r V Co aoco / 



For a related discussion see |67]- 
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We can now evaluate (37) for the Lifshitz-like solution, for which the relevant near-horizon 
data is given by 



2z 



^0 



l(z + 2)(z + 3)^^ 



„2 + 3 



and we recall that $ 



(z-l)(z + 3), 
and z = 



Co = 



24 



(z + 2) 



1 , /L6$ 

— m 
2a 



3 

2^' 



(38) 



6L-*(2-l) 



the scaling solution becomes: 



7? 



1 

47r 



24(1 



371 ' 



z + 2 



he 



We find that the viscosity to entropy ratio for 



7l<Ph 



+ 



12(z- 1) 
z + 2 



(326 



(y2-2a)(l>h 



(39) 



Let's analyze this result for the simplest case where we fix the dilaton couplings such that 
the global symmetry Q is preserved, i.e. 71 = and 72 = 2a. In this case (39) reduces to 



77 
s 



1 

47r 



(40) 



z + 2' z + 2 

Thus, for this special choice of 71 and 72 the correction to r]/ s takes on a particularly simple 
form^^ dictated by the value of the dynamical exponent z, which we recall is determined by the 



dilatonic coupling a. Note that when z = 1, (40) reduces to the well-known case [^3 EOj of an 
uncharged AdS brane, 

1 



s 



Air 



1 



(41) 



In that setup, the coupling f3i of the higher derivative terms could be related to the central 
charges a and c of the dual 4D CFT by making use of the trace anomaly^'', which gave /3i oc 



Thus, (41 ) could be expressed entirely in terms of the parameters of the UV CFT. For the present 



case, in (40) we see that ry/s is no longer uniquely determined by the UV central charges (through 
its dependence on (3i) , but it is also sensitive to z, which parametrizes IR data. Thanks to the 



In the absence of higher derivative corrections, the result f = 3^ had already been obtained for systems with 
Lifshitz scaling at extremality in |68) . 

As a consistency check, we note that the /3i term here agrees with the result of 48 , which included only 
Gauss-Bonnet corrections and did not couple them to the scalar. 

Similar arguments using the R-current anomaly relate /32 to a and c in a very similar way. 
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different scalings of the solutions in our model, the expression for rj/s in (40) is explicitly sensitive 
to both the IR and the UV data. We will come back to this point in the next section. 

Another interesting feature becomes apparent when one notices that the dependence on z can 



be removed by making use of the expression (38) for the scalar field at the horizon^^. Trading 



the explicit z dependence for {Q, t/^, e"'^'*}, we find that (40) becomes: 



s 



1 



1 - 4 ( 2/3i - + 1/32) 



2 T2„-2a4 



Q'L'e 



3r^ 



(42) 



With the charge rescaled appropriately, Q = ^^'^ ^ ^ , the expression (|42|) becomes identical 



in structure^^ to the standard result^'' for a RN-AdS brane (with charge Q) in a theory with 
higher derivative corrections of the form 5C = jSiR^yp^jW^^"^ + I32R^upctF^'^ F^"' + . . ., originally 
found in |37| I56|. This is not too surprising, however, if we recall that 77/5 is a horizon quantity, 
and at the horizon the role of the dilatonic gauge kinetic coupling e^"*^ is essentially to rescale 
the charge. We emphasize that this is a consequence of having chosen the dilatonic couplings 
in such a way to preserve the shift symmetry. Thus, even though the low-T Lifshitz solution 
of our dilatonic system is entirely different from the RN-AdS brane of Einstein-Maxwell theory, 
the structure of rj/s - as long as the shift symmetry is not broken - will be the same. 

Finally, we would like to comment briefly on ramifications for the KSS bound. Assuming the 
/3i coefficients are strictly positive, we see that for the case where the dilaton coupling respects 
the shift symmetry, as in (40), the f3i term always decreases the value r]/s whereas the /S2 



term always increases it. Moreover, as we take z to be very large the /3i contribution becomes 
vanishingly small while the /32 term approaches a constant value. Thus, in the large z limit, 
the R^upijF^'^FP"' term dominates the correction. The same behavior was also observed in the 
RN-AdS black brane correction when evaluated at extremality \37\ l56] . 



Using e^""^'- 



and recalling the form of $, one finds 



1 + 



^2^2g-2o.0„ 



Note that in order to have z > we need < Gr^e^"'*'' /L^. 

We stress however, that even though the form of r^/s in (42 1 is suggestive, the charge Q of the RN-AdS brane 



would not be the same as the charge Q measured at infinity from the solution with the near-horizon Lifshitz 
scaling. 

Note that our gauge field normalization differs by a factor of 1/4 from that of [371 156) . 
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We close this section by commenting briefly on the case where the shift symmetry is exphcitly 



broken, and both 71 and 72 are ahowed to take arbitrary values. By examining (39) we see 
that turning on 71 gives additional freedom to change the sign of the correction - and the /3i 
contribution can naively be forced to vanish by appropriately tuning 71. Clearly the size and 
sign of the correction to the universal 1 /47r result will now depend on the ranges (and signs) of 
all the parameters {7j,/3j} controlling the higher derivative terms, as well as on the value of the 
dynamical exponent. Thus, once the shift symmetry is broken, the viscosity to entropy ratio 



extracted from the scaling solution (52) is no longer related in any simple way to the value it 
would take for an AdS brane. As we will see in the following section, this is tied to the fact that 
in this model t//s not only runs as a function of temperature from its IR value to that in the 
UV, but can also have additional dependence on T within the T <C /x regime itself. 

IV. TEMPERATURE FLOW: FROM LIFSHITZ TO ADS 

We are now ready to discuss the behavior oirj/s as a function of temperature, and in particular 
in the two opposite regimes T <^ n and T » /i, the IR and the UV respectively. To get a better 
feel for the structure of 7]/s, let's start by assuming that the scalar field does not couple directly 
to the higher derivative corrections, 7^ = 0, and turn off the R^upuF^'^F''"' term by setting 
(32 = 0, so that the global symmetry (Q of the action is restored. We then see from (37) that 
rj/s can be written in the simple form 



s An 



1 



(43) 



Co 

with the 0(/3i) higher derivative correction parametrized entirely in terms of cq - different near- 
horizon geometries described by different values of cq. Moreover, since the near- horizon behavior 



of the solutions in our model changes as T and /i are varied, (43) shows in a clean way that the 
viscosity to entropy ratio will flow as a function of T/fj, - the temperature dependence encoded 
in Co translating into a flow for rj/s. Let's describe this simple flow in a bit more detail. 

At low temperatures T <^ fi the near-horizon geometry probed by the shear fluctuation is 



that of the Lifshitz-like brane. Reading off cq from (38), equation (43) tells us that in the IR 
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the viscosity to entropy ratio has the remarkably simple form 



1 

47r 



1 



24 



■A 



(44) 



IR 47r [ 2 + 2 

which we could have equivalently expressed^^ in terms of the dilatonic coupling a. On the 
other hand, for temperatures much higher than the chemical potential, T S> /i, the near-horizon 
geometry becomes the AdS-black brane - the high temperature has essentially 'washed-out' any 
footprints of the Lifshitz region in the IR. In particular, when the temperature is so high that 
/i is negligible, the near- horizon behavior is that of an uncharged black brane in AdS, for which 



both the gauge field and the scalar are constants. For this geometry cq = i^/8, and (43) reduces 
to 



uv 



1 

47r 



[1 



(45) 



in agreement with the well-known result of |691 [70] . As a simple consistency check, note that ( 44 ) 



reduces to the uncharged AdS brane result ( 45 ) when the dynamical exponent becomes z = 1 , as 



it should. Finally, when the temperature is large compared to /x, but not so large that /x should 
be neglected, the geometry starts being sensitive to the presence of charge. Unfortunately, due 
to the non- vanishing scalar field in our system, the usual RN-AdS black brane is not a solution 
of the theory, and in this intermediate temperature regime the near-horizon behavior is not 
known analytically. As a result, r//s must be computed numerically. Understanding in detail 
this middle region is, however, beyond the scope of this work. 

Although we don't have an analytic expression for the near horizon geometry for arbitrary 
temperatures, we emphasize that the parameter cq varies smoothly (and monotonically) |51| 
between the value it takes in the IR and that in the UV, 



I? I? 

— < Co < — (z + 2). 



implying a smooth interpolation for ry/s as a function of temperature^^ : 



47r 



s 47r 



24 

1 /3i 

2 + 2 ^ 



(46) 



(47) 



■ft 



Trading z for the coupling a we have 7 = 5 
Note that this behavior is similar in spirit to the temperature flow seen in higher derivative theories at finite 
chemical potential (37| i56j . 
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As we mentioned in Section III the holographic Weyl anomaly relates the coupling of the 
curvature-squared correction ~ PiR^yp^R^^^^ to the central charges a, c of the dual 4D CFT: 



^1 oc 



(48) 



Thus, comparing the left- and right-hand sides of (47) we see that, while at high temperatures 



r]/s can be expressed entirely in terms of the UV central charges {a, c}, for low temperatures 
the correction is 'dressed' by the dynamical exponent z, and is therefore also sensitive to IR 
data. This feature sets this construction apart from the setups of e.g. [691 [70], where UV data 
was enough to specify r]/ s. A similar behavior for transport - with a sharply different rj/ s in 
the IR and in the UV - was also found in [S3] , where it was intimately tied to the presence of a 
superfluid phase transition. 



Given the relation (48), it is natural to try to describe the flow (47) at arbitrary temperatures 



in terms of some combination Cg// of effective central charges, encoding information about the 
number of degrees of freedom of the dual system: 



s 



1 
47r 



1-C, 



(49) 



Clearly, this function would have to reduce to our results at the endpoints of the flow, 



eff 



OC < 



(50) 



24 



where we can associate Pi with the appropriate central charges of the dual UV CFT. Notice 



Pi 



that, as long as z > 1, the deep IR value (44) oirj/s is always larger^^ than the result (45 ) in the 
UV, with a monotonic flow between the two. Thus, at least in this simple case, the size of Cg// 
decreases with temperature from the UV to the IR, suggesting that it may in fact be a measure 
of the degrees of freedom in the system. It would be interesting to make this type of relation 
more precise. 

Let's now go back to discussing the generic higher derivative action ([s]), with all the couplings 
{Piyji} turned on. First we note that for T ^ fi, corresponding to the case where the near- 



It is largest in the z — >■ cxa limit, in which it approaches l/4n. 
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horizon is AdS-like, the correction to t//s is only sHghtly changed from that in (45): 



1 

UV 47r - 



1 



he 



,7i</'0 



(51) 



where 4>o is the boundary value of the (constant) dilaton, as given in (13). We therefore focus 
the rest of the discussion on the more interesting case of the low temperature Lifshitz-like brane. 
For convenience, we repeat here the full expression for r//s in the T <C )U regime, 



1 

IR An 



1 



24(1 



z + 2 



(52) 



with arbitrary couplings, breaking the shift symmetry of the two-derivative action. Compared 



to what we had in (44), note that r/fs is now sensitive to the horizon value <ph of the scalar. In 



the Lifshitz-like regime, the gauge kinetic coupling at the horizon scales with temperature as 



Since /i is the only other scale in the theory, simple dimensional analysis tells us that 



(53) 



(54) 



is the appropriate scaling of the dilatonic factor with T and ^. For the precise proportionality 
constant, we would need the exact dependence of /i on the parameters of the near-horizon IR 
solution, which in our model is not known analytically |501 151). 



However, from the scaling (54) we can see that the dilatonic factors e^^'^'^ and e"^"'^'^ in ^ 



will lead to additional temperature dependence in the IR, and give rise to a non-trivial behavior 
for r]/ s. Even in the absence of direct couplings of the scalar to the higher derivative terms 
(i.e. 7.j = 0), and setting /3i = for simplicity, we already see a new temperature-dependent 
contribution: 



1 

IR 47r 



oc /32 



6/2 



(55) 



The simple /32 = limit of (52) also leads to additional temperature dependence, but is a bit 



more subtle. In this case the deviation of rj/ s from l/47r becomes: 

^ /2-\-^V3(^-l)/2 



IR AtT 



OC Pi 



(56) 
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For positive 71 this result should not be trusted at very small temperatures - near extremality 
the correction stops being perturbatively small, signaling a break-down of our analysis. This be- 
havior can be traced back to the dependence of the scalar field at the horizon on the temperature 
of the system. In fact, divergences in the extremal limit of Lifshitz solutions are well-known (see 
|7H [72] for a recent discussion of stability issues of Lifshitz backgrounds), and higher derivative 
corrections are expected to become large in the zero temperature limit. A more quantitative 



analysis of the break-down of ( 56 ) would require a careful study of the solutions in the regime 
of nearly zero temperature, as well as of the backreaction of the higher derivative terms. 

Here we will content ourselves with pointing out that the presence of generic (ph dependence 
in r]/s leads to an additional sensitivity to temperature, even in the low-T regime, which one 
would not see if the global symmetry was preserved. In fact, contrast the cases we have just 
discussed to the situation with 71 = and 72 = 2a. Now the shift symmetry is intact, any 



dependence of (52) on (ph is lost and the IR expression for r]/s reduces to 



with no trace left of the additional running with T. In summary, much like in [55], it is the 
breaking of a symmetry (a shift symmetry here) which plays a crucial role in generating a 
particularly non-trivial structure for the viscosity to entropy ratio. 

Sample plots of the temperature dependence^^ of 7]/s in the IR can be found in Figures [l] 
and[2j In Figure [T] we have chosen to display the simplest dependence oi rj/s on T//i, in which 
all parameters except for z and /32 are set to zero. This case corresponds to that described in 



(55), and the resulting curves show a simple power law temperature behavior for {rj/s — l/Air), 
with the power determined by the value of the dynamical exponent z. Thus, one can control 
how sharply r]/s varies with temperature by appropriately tuning z. 

Going back to the more general case described by ( 52 ) , in Figure [2] we have fixed instead 



z = 3 as well as the values of /3i and /32, but have allowed 71 and 72 to vary between each 
curve. The parameters have been chosen to allow for minima for ry/s at some temperature, 



In the plots, we have absorbed the (unknown) constant of proportionahty appearing in (54 1 into the values of 
Pi and 132, since the latter are free parameters. However, we have insisted everywhere in the plots that /3i and 
132 are small. 
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C-l >5 0.082 



0.078 




FIG. 1: T/ jjL dependence of 77/s for a simple case with = 71 = 72 = 0} and /32 = 0.01. The plots are 
for values of 2; = {12, 9, 6, 3, 1}, from top to bottom. 

as observed in many substances in nature. We emphasize that while the presence of extrema 
is a rather typical feature^^ of the T/// dependence of our expression (52) for 77/5, the plots 



we have shown probe only a very small region of parameter space, and that a much wider 
range of behavior is possible. Clearly, for realistic phenomenological applications it would be 
valuable to construct more concrete models in which at least some of the parameters governing 
the temperature dependence can be fixed unambiguously. This is an important question that 
we plan to address in future work. 

V. DISCUSSION AND SUMMARY OF RESULTS 



A question that has emerged quite naturally in the context of the holographic gauge/gravity 
duality is whether hydrodynamic transport coefficients depend on the radial direction of the 
gravitational background used to compute them. On the field theory side of the duality, the 



^® Extrema can be easily obtained by choosing appropriately the signs and exponents of the competing contribu- 



tions of the j3i and /32 terms in ( 52 1 
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FIG. 2: T / IX dependence of rj/s illustrating existence of minima, reminiscent of data for various liquids 
in nature. All curves have the same parameter values of z = 1, /3i = 0.01, and /32 =0.1 with differing 
values of 71 = {-2.5, -2, -1.5, -1} and 72 = {-0.2, -0.2, -0.25, -0.75}, from top to bottom. 

radial dependence is understood as dependence on the RG scale. It has been known for some 
time that - in the AdS/CFT framework - the shear viscosity 77 can be extracted from quantities 
which do not depend on the holographic radial coordinate. In other words, it does not run with 
energy (see [391 El] for related recent discussions) . 

Even though rj/s doesn't exhibit RG flow, its behavior is not entirely trivial. It is possible 
to have constructions in which rj/ s takes on very different values at low and high temperatures 
(which we refer to as the IR and UV, respectively) , with a non-trivial temperature flow connect- 
ing the two regimes. One way of achieving such a flow in rotationally symmetric backgrounds is 
to include higher derivative corrections in the theory, coupled to a scalar field. This was done for 
example in p5], where the temperature dependence of ry/s - and the fact that its IR behavior 
was very different from that in the UV - was linked to the presence of a non-trivial scalar field 
controlling the strength of the higher derivative operators, as well as to the presence of a phase 
transition. It is natural to ask, then, whether the running of r}/s with temperature, and the 
type of "decoupling" of IR from UV physics observed in ^55j, can be seen in other setups, and 
not necessarily involving phase transitions. 
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With these motivations in mind, in this paper we have chosen to work with a model of 
Einstein-Maxwell-dilaton gravity, whose black brane solutions interpolate smoothly between a 
near-horizon Lifshitz scaling in the low-temperature regime and standard AdS scaling at high 
temperature. The temperature dependence of the solutions translates into a temperature flow 
for r//s, once higher derivative corrections are turned on. Thus, this construction provides an 
explicit setting in which r//s exhibits a non-trivial behavior as it varies as a function of T//i 
from the IR to the UV. To illustrate some of the features with a simple example, we summarize 
the case of /32 = 71 = 72 = 0, for which: 



1 
s 



^[l- 8/3i] , T » /i (UV) , 

(58) 



47r 



1 24 o 



T < /X (IR) . 

When the chemical potential is the dominant scale in the system, the putative field theory is in 
the non-relativistic scaling regime, hence the IR dependence of r]/s on the dynamical exponent 
z (or equivalently a). In the relativistic high temperature regime, on the other hand, any 
knowledge of z is lost - the shear fluctuation is now probing a geometry whose near-horizon 



behavior is AdS-like. The qualitatively different behavior seen in (58), depending on the relative 
hierarchy between the two scales in the theory, is reminiscent to that of [55] . 

The simple high-T result rj/s = 3^ [1 — 8/3i] matches that obtained from curvature-squared 
corrections in AdS [69j. In that case, the couplings of the higher derivative terms could be 
related to the central charges of the dual UV CFT via the holographic Weyl anomaly (in our 
notation /3i = |^^), so that one could write rj/s entirely in terms of UV data. The situation 
in our model is quite different: although at high temperatures it is analogous to that of [69j, in 
the low-T regime the behavior of r//s is no longer determined only by UV quantities, but it also 
depends on IR data - in this case, the dynamical exponent z. 



In addition to the 'simple' T/fi flow that is visible from (58), describing how r//s varies 
between its low-T and high-T values, there is another source of temperature dependence, inti- 
mately tied to the presence of the running scalar field. In fact, in the Liftshitz scaling regime 
simple dimensional analysis arguments tell us that 

e-^"'^" ~ i-j , (59) 
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with the horizon value of the scalar. Thus, any factors of (ph in the viscosity to entropy 
ratio will lead generically to additional sensitivity to temperature. We can see explicitly the 
appearance of these new temperature terms by looking at the form of r//s in the T <C ^ regime, 
for arbitrary higher derivative couplings: 
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IR Ait 



1 



24(1 - 

z + 2 



1 "ll 7 



1) 



z + 2 



/32e 



(72-2o)(/>h 



(60) 



The factors of will source new temperature dependence through ( 59 ) , even when the scalars 
are not directly coupled to the higher derivative operators (7^ = 0). 

Interestingly, the terms that generate this additional T//i dependence disappear when the 
higher derivative couplings are chosen to satisfy the shift-symmetry 



+ 5, 



A 



'A^ 



(61) 



of the two-derivative action. Specifically, when the dilatonic couplings are such that the sym- 



metry is preserved, i.e. 71 = and 72 = 2a, the expression (|60|) reduces to 

24 
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IR AtT 
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z + 2 



^ z + 2 ^ 



(62) 



with any trace of the additional temperature factors lost. Moreover, as discussed in detail in 



Section III A, the structure of (62) turns out to be equivalent - after an appropriate rescaling 



of the charge - to that of a RN-AdS brane. While this may seem surprising at first, it is easily 
explained by noting that when the shift symmetry is preserved, the Lagrangian takes the simpler 
form 

C = R-2A-2 {VcPf - e^^'t'F^ + L^/^i R^up.R^"^" + L^P2 e^'"^>Rf,^.p^F^"' PP" . (63) 

Thus, the only role of the dilatonic coupling e'^°"^'^ at the horizon - where rj/s is evaluated - is 
to rescale the charge of the U{1) gauge field. 

This construction has allowed us to explore some of the features needed to generate a tem- 
perature flow for the transport coefficients of a strongly coupled plasma. The key ingredients in 
our model were a non-trivial scalar field profile and different scalings in the IR and in the UV. 
The dynamical exponent z and the running of the scalar field conspired to trigger temperature 
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flow, and to generate a particularly non-trivial structure for 77/s, with the breaking of the shift 
symmetry playing a crucial role in giving rise to an additional source of temperature depen- 
dence. Although it would have been interesting to understand the intermediate temperature 
regime, here wc have restricted our analysis to the cndpoints of the flow, since this was enough 
to make the points we were after. Finally, we should note that it would be valuable to extend 
this work by performing a more systematic - and more quantitative - study of the temperature 
dependence of the transport coefficients of a strongly coupled plasma, especially in view of the 
current applications to the physics of the strongly coupled QGP. With this context in mind, 
we emphasize that the key ingredient for achieving non-trivial temperature dependence in the 
IR was the existence of a running scalar, and not necessarily the presence of non-relativistic 
scaling. Thus, similar constructions should conceivably achieve an analogous temperature flow 
in relativistic settings more applicable to QGP-like plasmas. 
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